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Abstract 

The finite g-oscillator is a model that obeys the dynamics of the 
harmonic oscillator, with the operators of position, momentum and 
Hamiltonian being functions of elements of the (/-algebra su 9 (2). The 
spectrum of position in this discrete system, in a fixed representation 
j, consists of 2j+l 'sensor'-points x s = ^[2s] q , s G {—j, ■ ■ ■ ,j}, 

and similarly for the momentum observable. The spectrum of energies 
is finite and equally spaced, so the system supports coherent states. 
The wave functions involve dual g-Kravchuk polynomials, which are 
solutions to a finite-difference Schrodinger equation. Time evolution 
(times a phase) defines the fractional Fourier-g-Kravchuk transform. 
In the classical limit q — > 1 we recover the finite oscillator Lie algebra, 
the N = 2j — > oo limit returns the Macfarlane-Biedenharn g-oscillator 
and both limits contract the generators to the standard quantum- 
mechanical harmonic oscillator. 
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1. Introduction 

Discrete models which are counterparts to well-known continuous sys- 
tems, and in particular those which contract to the standard harmonic os- 
cillator, are of fundamental interest in theoretical physics [l]-[5]. Moreover, 
finite discrete models are also of interest for the parallel processing of sig- 
nals through nano-optical devices, where the input may involve lasing carbon 
tubules, the output being registered by a finite sensor array, and the device 
consisting of a shallow planar waveguide — an oscillator which can carry only 
a finite number of states [6]. The salient purpose of such a device is to 
perform a finite analogue of the fractional Fourier transform [7]. 

In previous works on the one-dimensional finite oscillator [3, 6, 8], oscilla- 
tor systems were characterized in the familiar context of Hilbert spaces and 
Lie algebraic theory; in [5] these requirements were formalized into the three 
following postulates: 

1. There exists an essentially self- adjoint position operator, indicated Q, 
whose spectrum T,(Q) is the set of positions of the system. 

2. There exists a self-adjoint and compact Hamiltonian operator, H, which 
generates time evolution through the Newton-Lie, or equivalent Hamilton-Lie 
equations: 

iH,iH,Q]]=Q ~ : : -£ (i) 

where [ • , • ] is the commutator. The first Hamilton equation in (1) defines the 
momentum operator P, while the second one contains the harmonic oscillator 
dynamics. The set of momentum values of the system is the spectrum E(P) 
of P. 

3. The three operators, Q, P and H, close into an associative algebra, i.e., 
satisfy the Jacobi identity, 

[P, [H, Q]\ + [Q, [P, H]] + [H, [Q, P}} = 0. (2) 

The second and third postulates determine that [Q, P] must commute 
with H, which implies that it can only be of the form [Q,P] = iF(H), 
where F is some function of H (including constants) and the i is placed to 
make F(H) self-adjoint, but do not otherwise specify this basic commutator 
further. For a constant F(H) = hi, one recovers the standard oscillator 
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algebra H4 = span {H, Q, P, 1}, which contains the basic Heisenberg-Weyl 
subalgebra W\ = span {Q, P, 1} of quantum mechanics. In our first works 
[6, 8] we examined the cases which, in the unitary irreducible representations 
of spin j = t;N (N G {0,1,...} fixed), correspond to the linear function 
F(H) — H — {j + |)1 =: J 3 , and so the operators close into the Lie algebra 
so(3) = su(2) = span {Q, P, J 3 }. The purpose of the present paper is to 
study the case when, for q := e~ K , the basic commutator is 

[Q,P] = iF q (H), H=J 3 + (j + ±)l, (3) 

F(H) = c~ 2 " J * C ° Sh ^ c -"-/3 COsh (j + ^)« U) 
qK ' 2sinh|/t 2sinh|/t v ; 

= \ e - Kji (e- Rjs cosh \n-T 2j+l (cosh / sinh \k, (5) 

where T n is the Chebyshev polynomial of the first kind, and q G (0, 1] or 
k G [0, 00). In particular, F±(H) = J 3 returns the previous su(2) case [6]. 

An important ingredient for the postulates of harmonic oscillator dynam- 
ics is an unambiguous correspondence between the physical observables of 
position, momentum and energy, with the elements of the associative alge- 
bra. In Section 2 we recall the main relevant results on the algebra su g (2) and 
its standard representation basis. The su g (2) nonstandard basis, investigated 
in [9, 10], is introduced in Section 3 to exhibit our proposed correspondence 
explicitly in terms of the generators of su ? (2). With our postulated choice, 
the position and energy spectra in the (2j'+l)-dimensional representation 
j = \N of su g (2) will be 

S(Q) = x s = \[2s] q = 2 S ^T ft , { ./•-./• 1 ./} -' V W 

E(if) = E n = n+\, ne{0,l,...,2j}=:n\ 2 J . (7) 

We recall the definition of the q- number for q = e~ K : 

q\ r — q~\ r sinh \tk , , 

\r\ q = [71,-1 = -l-r] := — = . \ ■ 8 

q y H qi — q~2 smh |/c 

Note that the g-number of an integer r is C/ r _i(cosh the Chebyshev 
polynomial of the second kind. The spectrum of momentum is the same 
as that of position, S(P) = T,(Q). The classical limit is lim^i [s] q = s, 
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when the g-algebra su q (2) becomes the Lie algebra su(2); then, the set of 
positions become equally spaced and we are back at the previously known 
finite oscillator [6]. But for all other values of the deformation parameter 
q, the 'sensor points' of the system are concentrated towards the center of 
the interval, while the endpoints are spread farther apart. Yet the energy 
spectrum remains an equally-spaced set, and therefore the system follows 
harmonic motion. 

The finite g-oscillator wave functions are the overlaps between the position 
and energy eigenbases. They are written out in Section 4 in terms of the dual 
g-Kravchuk polynomials, and are orthonormal and complete over the sensor 
points of the system. The momentum representation of these wave functions 
is addressed in Section 5 with the Fourier-g-Kravchuk transform, and in 
Section 6 this transform is fractionalized. The evolution in time of a finite 
g-oscillator (or equivalently, the parallel processing of a finite signal along 
the axis of a shallow planar waveguide), is the 2- fold cover of the fractional 
Fourier-g-Kravchuk transform matrix; the metaplectic sign appears thus for 
half-integer values of j, which corresponds to a finite systems of an even 
number of points. In Section 7 we introduce the concept of an equivalent 
potential for discrete systems which is based, as in the continuous case, on 
the existence of a ground state with no zeros. Finally, in Section 8 we verify 
that the contraction limits q — > 1 and N = 2j — > oo of the algebra su g (2) 
reproduce the known results for the finite oscillator and the continuous q- 
oscillator. The corresponding limits for the wave functions however, present 
further challenge. 

2. The algebra su g (2) and its standard basis 

The quantum algebra su q (2) is the associative algebra generated by three 
elements, usually denoted as J+, J_, J 3 , subject to the commutation relations 

[Js,J±] = ±J±, [J + ,J-] = [2J 3 ] q . (9) 

Equivalently, writing J± — J\ ± iJ 2 , we characterize the algebra su g (2) by 

[J2,Js]=iJi, [J 3 ,Ji]=iJ 2 , [Ji,J 2 ] = §[2J 3 ] ? . (10) 

The first two commutators in (10) have the structure of the oscillator Hamil- 
ton equations (1), while the third one involves the g-number (8), which dis- 
tinguishes g-algebras from Lie algebras, the latter corresponding to the case 
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q = 1. The following element in the covering algebra of su g (2) commutes 
with all others, 



C q := Jl + Jl + [h-\f q + \[2M 

112 _ 1 

2'q 4' 



= J + J_ + [J 3 - 112 - 1 



(11) 



and is called its Casimir operator. 

It is convenient to have a realization of the su g (2) generators in terms 
of first-degree differential operators, acting on spaces Tij of functions of a 
formal variable x, and depending on the numerical irreducible representation 
label j. This is 



J+ ■= Ji + iJ 2 <-> x[2j-x-^-] =x[j-J 3 ] q , 



dx 



J- := Ji — i J2 



ax 



— \ X ~T~] =-[3+Mqi 

x l dx 1 i x 
3 E {0, §, 1, . . .} fixed. 



(12) 
(13) 
(14) 



The set of power monomials x^ +m \ 3 m= _j are eigenfunctions of J3 and pro- 
vide the standard basis for the irreducible space Tlj, of finite dimension 2j+l. 
The functions of the basis were chosen in [9, 10] with the following constants: 



P m {x) := q* 



2j 
j+m 



1/2 



x- 



.j+m 



(15) 



where the g-binomial coefficient m is defined (using the standard notation 

L n J q 

of g-analysis [11]) for m > n nonnegative integers by 



rn 
n 



n-l 



-i)V 



(^9)n := n(!-^), n = l,2,3,..., 



i-|n(n-l) (g ""j?)" 
(95 9)n 

9)0 = 1. 



(16) 
(17) 



fc=0 



For any two complex vectors a, b e H,, 



^ a mf m , b — ^ ftmfm, 



m=-j 



m=-j 



(18) 
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there is a natural sesquilinear inner product 

(a,b) Wj := ]T a* m (3 m , (19) 

m=-j 

with respect to which the standard basis is orthonormal. The action of the 
su g (2) generators and Casimir operator on the standard basis is well known: 

■hfL = mfL J±fL = J\j±m+l] q \jTm] q f m±1 , (20) 

Cqf m = C q f 3 m) C q \= [j + ^lq ~ 4- (21) 

These equations are of course independent of the realization of the basis 
vectors by the power monomials f m (x) in x. 

The spectrum of the diagonal generator J 3 [see (14) and (20)] is linear 
and bounded, as that of a finite version of the quantum harmonic oscillator. 
Indeed, this is our choice for the finite g-oscillator Hamiltonian, displaced so 
that the ground state has energy |, namely 

H = J 3+3 + l Hfi = (n+l)fi, n:=j+m, (22) 

where is the mode number that counts the number of energy quanta. 
At this point we are presented with what would appear as a 'natural' as- 
signment for the position and momentum operators, Q «-> J\ and P <-> — J 2 , 
because it would be the simplest generalization of the previously studied 
q = 1 case [6, 8]. This choice would bring the first two commutators in 
(10) to reproduce correctly the two Hamilton equations in (1), while the 
third commutator [Q,P] would have the form (3) with F q (H) = |[2J 3 ] ? = 
sinh k(H — j — |)/2 sinh \k. In this 'naive' model however, the spectra of Q 
and P are not algebraic; they must be computed numerically as roots of a 
polynomial equation of degree 2j + 1. 

3. The nonstandard basis 

While we do not discard the model suggested at the end of the previous 
Section, we find more attractive to propose a correspondence between the 
physical observables of position and momentum, Q, P, and the nonstandard 
(also called twisted) operators (see [12]-[14], [9, 10]), which have the virtue of 
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possessing an algebraic spectrum x s := |[2s] g , s\ 3 _j- The position (and hence 
momentum) observables will be thus identified with the following operators: 

Q = J, ;=gi J 3 Jig iJ 3; (23) 
-P = J 2 :=qi J *J 2 q^, (24) 

while the Hamiltonian H is associated to J 3 by (22) as before. 

We note that while the g-number (8) displays symmetry under (/-inversions, 
q <-> q^ 1 , [r] q = [r] ? -i, the identification of tilded operators in (23)-(24) pre- 
serves this symmetry with the concomitant reflection J 3 <-> — J 3 . This means 
that the ground state of a q < 1 oscillator is the top state of its q^ 1 > 1 
partner. 

The commutation relations among the nonstandard operators and J 3 are 

[J 3 ,Q] = -iP, [J 3 ,P] = iQ, (25) 
[Q,P] = ^q^{q-hj + J_-q\j_J + )q^=:\F q {C q ,J 3 ) (26) 
= i(e- K - h [(C q + \) sinh i« + |csch \k] - \ e~ 2Kj3 coth i/c), 

where g = e~ K as before. The operator F q (C q , J 3 ) defined in (26) commutes 
with J 3 and is also diagonal in the standard basis; in the irreducible repre- 
sentation j, 

. _ e" 2mK cosh|/t - e" mK cosh(j+i)/t . 
qJm ~ 2sinhi K /m ' 1 7j 

but its spectrum is not a good candidate for an oscillator Hamiltonian, be- 
cause it is not equally spaced [unlike (7)], and so the motion would not be 
harmonic, but dispersive. In terms of the position and momentum generators 
(23)-(24), the Casimir operator (11) acquires the form 

C q = sechlK(Q 2 + P 2 )e Kj3 + D q (J 3 ), (28) 
D q ( J 3 ) := sech \k ([ J 3 -|] 2 - \e~ Rj3 coth i/c + ^csch \k) - \. (29) 

We recall a previous phase-space picture for the finite oscillator of 2j+l 
points, considered in [15], as the (classical) sphere Q 2 + P 2 + Jf = j(j+l), 
having circular sections of square radius Q 2 + P 2 = (j+|) 2 — (J3— \) 2 — ^3- 
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For su g (2), the corresponding surface now has the section 

q 2 +p 2 = ([j+w q cosh^-[j 3 -ir q 

+ \e- Kj ' A coth \k-\ csch \n) e~ Rj \ (30) 



Phase space for the finite g-oscillator is suggested thus as g-dependent pear- 
shaped spheroids, tip-up for q < 1 and tip-down for q > 1 (recall the q <-> q^ 1 
symmetry with the inversion of J3). The g-harmonic oscillator evolution (i.e., 
a phase times the so-defined fractional g-Fourier-Kravchuk transform) will 
rotate this space around the J3 vertical symmetry axis of the spheroid. 

In this finite g-oscillator model we interpret the eigenvalues x s of Q := J\ 
as the discrete values of the position observable. The eigenfunctions g\(x) 
and eigenvalues of this nonstandard operator were found in [9] , and they are 
of the form 

Qgi(x) = x s g J s (x), x s = \[2s] q = ^^T^, = ~ x -s, (31) 
gifx) = 1 i{q^ l - 2]) x-q)^ s (-q^ l - 2]) x-q) j+s = gl s {-x), (32) 

(33) 



Ts ■= Q' 



2j 

j + s 



-2s 



1 + q 
q i2(-q;q) 



2j 



They are normalized with respect to the inner product (19), and are orthog- 
onal because they correspond to distinct eigenvalues x s . This basis of 2j + 1 
functions g 3 s (x), s\ 3 _j we call the position basis. A signal consisting of 2j + 1 
values $ s , sensed at the positions x s [given in (6)], is 

*= ES.rf'eW,-, (34) 

s=-j 

and can be realized either as a function of x, or as a (2j + l)-dimensional 
column vector with components numbered by s\ J _j. 

4. Finite g-oscillator mode wave functions 

We have now two bases for Ti.f the standard basis {/4}m=~i °f fnode 
n = j + m (and energy E n — n + ^), and the nonstandard basis {giY s =-j op- 
position x s = \[2s] q . In the realization of su g (2) generators given in (12)-(14), 
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the mode basis is realized by the power functions in (15), and the position 
basis by (32). We can use this realization to find the unitary transformation 
between these two orthonormal bases, and thus define the finite (/-oscillator 
wave functions by the overlap 



n 



2j 

o > 



of mode n = j + m. 
on points x s = \[2s] q , s\ J _j. 



(35) 



By construction, this set of functions is orthonormal and complete under the 
Tij inner product (19). 

The overlap (35) is obtained by expanding the function g{{x) of (32) into 
a power series in x, which is then 



9 3 s (x) 



E *%i\xsy f^), E QfSSMnLw- (36) 



mi) i 



m=-j 



s=-j 



The expansion of g 3 s (x) in x is [9] 



gi(x) = 'f a E q^ j+m)U+m - 1} 



m=-j 



2j 
j+m 



1/2 



K j+m {\(j-s)--l,2j\q) P m {x), 

(37) 



expressed in terms of the dual q-Kravchuk polynomials, 



^„(g" 5 W~ 2j ;c,2j|g):= 3 02 



' q n , q ? , cq^ 2j 



q; q 



(38) 



where 302 is the basic hypergeometric function defined in [11], and the coef- 
ficients 7^ are given in (33). 

In the particular case of our concern, the argument of the dual g-Kravchuk 
polynomial is A(£) = q~^ + cq^~ 2 ° with c = — 1 in (37), is given in terms of 
the positions x s = \[2s] q , s | - , of the finite g-oscillator by 



X(j-s) = q-i^-q-i- 



-2e JK sinh ks 
2q~ j ~^(q -l)x s — -(Ae jK sinh \k) x s , 



(39) 



and q = e K as before. From (37) thus, the finite g-oscillator wave functions 
of mode number n = j + m, n\^\ are 
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q \(j+s)+\n(n-l) 





■ 2j - 




2f 


l+q~ 2s 


\ 


-j + s. 


g 2 


. n _ 


g2( y -q;q) 2j 



xK n (2q-*-*{q-l)x 8 ; -l,2j\q). (40) 
The explicit expression for the dual g-Kravchuk polynomials in this case is 



K j+m (X(j-s)] -1, 2j | q) = 302 



q-i- m ,q s -i,-q-i- s 



q; q 



q-^, 

^ (g~ j ~ m ; q )k(q- j+s ; q)k(-q- j - s ; q)k q k 

k=0 



(41) 



(<r 2j ; ?)* 



(e;9)fc' 



where (z;q)k is defined in (17). 

The lowest mode of the oscillator is [see (40) for n = j + m 



0], 





■ 2j " 


1 + q~ 2s 


\ 


J + s. 


<? 2 2(-g; g) 2 j 



7s- 



The finite g-oscillator wave functions possess definite parity, 

aSW-s.) = = (-i) n $i 2j|9) (^), 



(42) 



(43) 



and, as is to be expected, in the limit q — > 1 return the Kravchuk functions 
of the finite oscillator [6] 



lim$i 2 ^(:r s ) = 2' 



2j 
j + s 



K n (j-s;±,2j), 



(44) 



with the classical Kravchuk polynomials, introduced by Kravchuk in [16]. 

The dual g-Kravchuk polynomials - as all polynomials - are analytic 
functions on the complex plane of their argument. As before in the finite 
oscillator models [6, 8, 17], this argument is the position coordinate, which 
can be analytically continued to real or complex values X, even if the inner 
product of the space Hj is only over the point set {x s }, s\j. As to the 
g-Kravchuk wave functions (40) the factor in front of the polynomial is a 
function that is analytic in the argument s within the interval — j — 1 < s < 
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j + 1; this means that in the position coordinate, analytic continuation is 
possible within the interval x-j-i < X < Xj + \. 

5. Fourier-g-Kravchuk transform to momentum space 

The identification of the position and momentum operators, Q — Ji, 
P = — J2 in (23)-(24), brings formulae (25) to the role of the two Hamilton 
equations (1). [This also holds for the 'first' choice using the standard basis, 
Q <-> Ji, P <-> — J 2 , that we outlined in Section 2, as well as for all oscillator 
models, finite or standard.] The evolution of the finite g-oscillator over time 
in quantum mechanics, or along the optical axis in the waveguide model, is 
thus the harmonic motion 

-\th(Q\wh . ( Q(T)\ ( cosr sinr\/Q\ . . 

6 \p) e = -{p { r } ) = {- S mr cosr ) [p] " (45) 

This is a U(l) group of inner automorphisms of the su g (2) algebra, and of 
rotations of the phase-space surface around its vertical axis. It covers twice 
the SO(2) cycle of fractional Fourier- q-Kravchuk transforms, /Q, of power 
a = 2t/tt and angle r, 

K* := exp(-i7ra(J 3 + j)/2) = e ina/4 exp(-i ixaH/2). (46) 

For a = 1 we have the Fourier-g-Kravchuk transform )C q . The action of 
Kq on the eigenbasis of position yields the eigenbasis of momentum, 

gi(x):=lC q gi(x). (47) 

These functions have the properties and form 

Pglix) = -Y r gUx), Y r = \[2r\ q = - ~Y- r , r\ j _ p (48) 

g 3 r {x) = g 3 r (\x) = gl r (-ix) (49) 

)j+ri 



7^'(i?4 (1 ^x; g)j_ r 2j) x; q). 



where 7^ is the constant given in (33); the spectrum of momenta, Y r , r\ 3 _j, is 
the same as that of position [cf. (31)]. Since Kft is unitary under the inner 
product in Tij, the Fourier-g-Kravchuk transform of the finite g-oscillator 
eigenfunctions (35)-(40) of modes n = j + m, are 

¥^\x s ) := K q ^\x s ) := (gi, JC q P m ) Hj = (-^(i,), (50) 
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as in all oscillator models. 

The Fourier-g-Kravchuk transform of a function or signal $, of values 
<&(x s ) = (g{, $)h- 011 the finite, discrete sensor point set {x s }, s|j, is defined 
by 

$(x r ) = S) Wj = E IfJfW $(x s ), (51) 

s=-j 

where the kernel is the overlap of the position eigenfunctions g J s in (32) with 
the momentum eigenfunctions g 3 r in (49), 

K(f q) :=(g,gi) Hj . (52) 
This kernel is given explicitly below in (56) with a = 1. 
6. Fractional Fourier-g-Kravchuk kernel 

The Fourier-g-Kravchuk transform (50) is fractionalized by the operator 
ICg in (46), independently of the realization, on the mode eigenbasis of J 3 , 

K? q fi = exp(-i7ra(j + m)/2) f' m = exp(-i7ma/2) P m . (53) 

When we apply K. q on a finite, complex 'signal' function of 2j + 1 points, 

*(x.) = (si, *)„ = ± (gi, fi) nj (fL *)« (54) 

m=-j 

we obtain another such function, labelled by a, 

& a \x s ) := /CJ*(s.) := {gm®) H] = (}C- a gi,^) nj 

= E (icrni,ti)H i (fL*)H i = E {gi^qfmMfL^n, (ss) 

m=—j m=~j 

= E e- i7ra(J+m)/2 (^',/4)^ E (/^rfWrf.*)*^ E k:^" 'Im.'v). 

m=~j s'=-j s'=—j 

where the fractional Fourier-g-Kravchuk transform kernel K^f J ^ is a (2j+l) x 
(2j+l) matrix of elements given by the bilinear generating function [18, for- 
mula (8.15)] 

tf^'l*) := E $i 2jl9) (^) e" 1 W2 4 2il5) (av)* (56) 

n=0 
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= J J sriM)sW 7 (-q- 2 J-H; tf-^-q-i-rf-^-q-i-^-t- q,-t), (57) 
where 

t := e~ [na/2 , (58) 

R (f] = (g a ~~ J *, -q- J -% q s '- J t, -g-i- s 't, -t; g)oo 

Ps ' S ' [ ) • (f-'t,-q*'W-%-q-'-'t,-q-Vt ] q) 00 ' 1 J 

and 7^ is given by (33). The function sW-i, defined in [11], is 

w , , , , v ^l-ag 2fc (a,b,c,d,e,f;g) k z k 

8 W 7 {a; b, c, rf, e, /; g, z) := — — rr^, 

^ I- a (q,qa/b,qa/c,qa/d,qa/e,qa/f;q)k 

(60) 

where (a, . . . , c; q)^ := (a; q)^ . . . (c; q)^ and (a;g)oo = U^oi 1 - °<? fc ) in 
accordance with (17). This function can be expressed in terms of the basic 
hypergeometric function §07 (see [11, §2.2, formula (2.5.1)]), with coefficients 
which allow it to be reduced to the basic hypergeometric function 4 3 : 

tWri-q-V-H; q s ~ j , -q- j ~ s , q s '~ j , -q'^ 3 ' , -f; q,-t) 

(-q- 2j t,q-l- a ',-q-t + '',t;q) 00 x ( q~ j+s ' , -q- j ~ s ' ,t, -t 



(-q-^% q-^'t, g-V, -t; q)^ \ -q^% q~^, -q «' ' ^ 

We also note that due to relation (a; q) n = (a; q)oo/(aq n ] q)oo, the expression 
for P s ,s'{t) in (59) can be reduced to 

R (A (q^t; q^-Aq-^'t; q^j-q-^t; g) J ,[ g ' »/: q) j+2s+8t 

(-9 3t ;qhj 

(62) 

Naturally, /Q 1 1C a q 2 = lC a q l+a2 and K? q = 1. The 'phase correction' by 
na = 2t which we introduced in (46) implies that K, q — 1 (as the ordinary 
Fourier integral transform), while the fourth power of the oscillator evolution 
operator exp(iriJ) is —1 for the full rotation angle r = 2rr. This is the 
analogue of the metaplectic sign of the waveguide case (see [19], cf. [7]), 
where the U(l) subgroup generated by the latter covers twice the SO (2) 
of the former. Parity is conserved under the fractional Fourier-Kravchuk 
transformation because J 3 commutes with the inversion of phase space. And 
again, in the limit q — > 1 we recover the previous Fourier-Kravchuk kernel 
expressed in terms of the Wigner 'little-cP functions [8], 

hmi^^ = = Hr *' d { s ,{\™). (63) 
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7. Equivalent potentials 



In ordinary quantum mechanics, the ground state ^o(x) of a system with 
a potential V(x) and energy Eq > — oo, has no zeros; thus, the Schrodinger 
equation determines the potential energy of the system from the ground state, 

(d 2 \ d 2 

-l-^ + V(x)-E \y (x) = V{x)-Eo = \—-V {x)/v Q {x). 

(64) 

As a well-known example we have the harmonic oscillator, whose ground 
state is fo(x) ~ e~^ x , so ^^0(2;) = (x 2 — l)^f {x) an d (64) yields correctly 
V(x) - E = \x 2 - \. 

In the case when the system is discrete over the set of points x s = sh + xo, 
with integer s, which are equidistant by h, an equivalent potential may be de- 
fined following (64). We qualify it as equivalent because the discrete systems, 
that have been studied (such as Kravchuk, Meixner and Hahn systems [l]-[4], 
[20, 21]), obey Schrodinger-type difference equations which do not separate 
into a sum of terms, where one is readily identifiable with the kinetic term 
of the second-degree difference operator, plus a potential term that is only 
dependent on position x s . The symmetric second-difference operator, acting 
on functions of be expressed in terms of the right-difference and the 

left-difference operators V R and V L , 



V R := = ^-(e* - 1) = -^—^ - 1), 

l\X s L\X S X s+ i — x s 

V L := ^ = -}-(l-e- d °) = —^(l- e - d °), 

\ X s \%s %s %s— 1 



(65) 



where A = e ds — 1 = e 9s V. So, a difference analogue of the differential 
operator d 2 /dx 2 in (64) has the form 

' -(V R -V L ). (66) 



Xs+l/2 — Xs-1/2 



Consequently, when the ground state of the system is ip(s) := ^/q(x s ), the 
equivalent potential, according to its quantum-mechanical correspondent in 
(64), is 

V(x s ) -E = —— T (V R - V L )^(s) (67) 

2^(S)[x s+ i/2 -Xs-1/2] 
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'•0(5+1) - 4>(S) l/)(S) — 1p(S — 1) 



2(^+1/2 -2^-1/2)^(3) 



Xs+l X s 



Xs ^5 — 1 



In the case of the finite Kravchuk oscillator, the set of values of position 
x s — s (h — 1 and xq = 0) is finite: {x s Ys=-j- Yet, the wave functions 
ip(s) := \I/q 2 (x s ) can be analytically continued in x everywhere except for 
branch-point zeros at x±(j + i) := ±(j+l), which are due to the square root of 
the binomial distribution. Thus, on the closed segment < x < Xj + ±, 

the second difference in (67) is defined for any real value of x in the interval 
X-j < x < xj. A similar extension and range of validity holds for the Meixner 
and Hahn oscillator cases [4, 20, 21]. The lowest mode of the Kravchuk 
oscillator, where h — 1, is given in (42). /,From this one derives the equivalent 
potential for the Kravchuk eigenfunction system 



V(x s ) -E + l 



if>(s + 1) + ip(s - 1) 
2ij){s) 

+ s)(j + s + l) + y/(j- s)(j - s + 1) 
2y/{j + 1) - s 2 



• (68) 



When the set of position values is not equally spaced, as is the case in 
the finite g-oscillator, {x s Y 8 =-j as in (31), we shall consider the differences 
with respect to the position coordinate 



X c 



sinh sk 
2l ~ ' 2sinh|/t 



x s +i — x s = cosh(s+i)fi;, 

x 



x s -i = cosh(s— ±)K. 



Taking into account that 



if;( s + l) = q ->>- 1 '\ 



cosh(s + 1)k sinh(j — s)k 



(70) 



\ coshs/t sinh(j + s + 1)k 

we arrive at the expression for the equivalent potential in the general case 

1 



V(x s )-E = 



2 q 1 / 2 cosh(s + \)k cosh(s — \)k 



x<q 



cosh(s 



cosh(s — 1)k sinh(j + s)k 
cosh sk; \ coshs/t sinh(j — s + 1)k 
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+ q 



s cosh(s — |)/c cosh(s + 1)k sinh(j — s)k 




(71) 



for functions ip(s) := W (x s ) (see formula (42)). Obviously, in the limit 
when q — > 1 (that is, k — > 0), this expression coincides with (68). 

Note that acceptable ground states occur for values of q which are larger 
of some number a < 1 (this number a changes with the value of j) while 
lower values of q present the raised-wings problem of interpretation. The 
corresponding potentials have an oscillator-type form for all values of q and 
this property is of course likewise shared by the g-Kravchuk wave functions. A 
study of these functions with attention to their oscillations and convergence 
should be undertaken but this task is beyond the purpose of the present 
paper. 

8. Contraction of the algebra su,j(2) — > osc q 

We consider a sequence of finite g-oscillators over sets of 2j + 1 points 
which increase in number and density as j — > oo, while the mode number 
n = j + m remains finite, i.e., near to the ground state n = (for eigenvalues 
m of J3 near to — j). The spectrum of the Hamiltonian operator H = J3+J + I 
of the g-oscillator retains its linear lower-bound spectrum (7) for all j's in the 
sequence. In the case of the (q = 1) finite oscillator, we showed in [22] that the 
su(2) dynamical algebra, wave functions, and Fourier- Kravchuk transform, 
contract to the ordinary oscillator algebra osc = span {Q, P, H, 1}. In the 
present g-case we follow an analogous contraction to the g-oscillator model of 
Macfarlane and Biedenharn [23, 24]; nevertheless, there are some important 
differences between the q- and non-g cases that we shall point out below. 

The 'sensor points' of our finite g-oscillator [i.e., the spectrum of Q e 
su g (2), S(Q) in (6)] extend between x_j and Xj, inside an interval which 
grows asymptotically with j as ~ q~i = e JK (for < q = e~ K < 1, k > 0) 
- and are not equally-spaced within. Our contraction process will keep the 
range of positions finite by introducing, for each finite j, the operators 



WjQ, P {j) :=WjP, 



(72) 
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scaled with coefficients whose asymptotic behavior is appropriate, 



Wj := - — ^— = e 2^+2> K . 



IX j 



\ H gV 2 (l - g) = e - JK Ve"2 K sinh ±k. 

\ smhj/t v v 2 

(73) 

The number operator, iV := H — \ = J 3 + j, is assumed to act on a subspace 
of functions whose mode eigenvalues n = j+m remain finite inn 6 {0, 1, . . .}. 

As we let j — > oo, the su g (2) algebra of the finite g-oscillator will contract 
to a different g-algebra, that will characterize the 'continuous' limit of our 
finite model. The commutation relations (25), which can be written 

[H, Q&] = -i P®, [H, P&\ = i Q®, (74) 

continue to be harmonic oscillator Hamilton equations. The third commuta- 
tor (26), which is characteristic of our su 9 (2) finite model, becomes 

[Q U \ P {j) ] = w][Q, P] = i q - F q (C q , J 3 ). (75) 

Xj 

Acting on the subspace of functions whose mode numbers remain finite, from 
(27) we find that the asymptotic behavior of the right-hand side of (75) is 

a<j+h) i 



Xj 



F q (C q ,J 3 ) ~ q J ^ = q H ~^ = q N . (76) 



When j — > oo, the formal limit operators — > Q and P^ — > P satisfy the 
oscillator Hamilton equations (74) and 

[Q,P}=iq N , N = H-\. (77) 

The reader may be more familiar with the contracted algebra span{Q, 
P, N} when it is written in terms of the raising and lowering operators as 

A± := Q T iP= lim J±, (78) 

whose commutation relations are 

[A + ,A_] = 2q N , A_A + -qA + A_ = l. (79) 
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This we identify as the g-oscillator algebra osc q defined by Macfarlane [23] 
and Biedenharn [24]. The j ; — > oo limit of (78) yields 

A + ^\X) = ^{n+l} q ^Ux), (80) 
A-^(X) = ^/W^SiPO. (81) 
where {n} q := (q n — l)/(q — 1) and 

¥«\X) = ^=(A + r¥ ( ?\x) (82) 
Vn! 

are mode eigenfunctions obtained from A_ty$\x) = 0. 

We would like to point out however, that before the limit j — > oo is 
achieved, the spectra of position and momenta, (31) and (48), are asymptot- 
ically constrained to a finite position interval 

|E(QW)| < WjX . „ 1/^/2(9-1 - 1). (83) 

Only in the q — 1 finite oscillator case [22], where = j, does the position 
interval grow to the real line as ~ yfj, keeping equal distances ~ 1/V7 
between neighboring sensor points. For any other < q < 1, all points a; s 
of E(Q < -' ? - ) ) except x± 3 -, will crowd towards zero in the middle of the interval. 
This feature of the contraction limit between g-algebras is at variance with 
that encountered with Lie algebras, where one can extend the operation 
from formal operators to finite Hilbert spaces of growing dimensions, to find 
limits from Kravchuk to Hermite functions, and Schrodinger difference to 
differential equations. This matter also requires a separate, deeper analysis 
that we leave for a separate publication. 

Acknowledgements 

We thank the support of the Direccion General de Asuntos del Personal 
Academico, Universidad Nacional Autonoma de Mexico (dgapa-unam) by 
the grant IN102603-3 Optica Matemdtica and sep-conacyt project 41051- 
F. A. U.K. acknowledges conacyt (Mexico) for a Catedra Patrimonial Nivel 
II. 



18 



References 

[1] Atakishiyev N M and Suslov S K 1989 A model of the harmonic oscillator 
on the lattice Proc. VI All-Union Colloquium "Contemporary Group 
Analysis: Methods and Applications" (Baku: Elm) 

[2] Atakishiyev N M and Suslov S K 1991 Difference analogs of the harmonic 
oscillator Teor. Math. Phys. 85 1055-62 

[3] Atakishiyev N M and Wolf K B 1994 Approximation on a finite set of 
points through Kravchuk functions Rev. Mex. Fis. 40 366-77 

[4] Atakishiyev N M, Jafarov E I, Nagiev Sh M and Wolf K B 1998 Meixner 
oscillators Rev. Mex. Fis. 44 235-44 

[5] Arik M, Atakishiyev N M and Wolf K B 1999 Quantum algebraic struc- 
tures compatible with the harmonic oscillator Newton equation J. Phys. 
A: Math. Gen. 32 L371-L376 

[6] Atakishiyev N M and Wolf K B 1997 Fractional Fourier- Kravchuk trans- 
form J. Opt. Soc. Am. A 14 1467-77 

[7] Ozaktas H M, Zalevsky Z and Alper Kutay M 2000 The Fractional 
Fourier Transform (Chichester: Wiley) 

[8] Atakishiyev N M, Vicent L E and Wolf KB 1999 Continuous discrete 
fractional Fourier transforms J. Comput. Appl. Math. 107 73-95 

[9] Atakishiyev N M and Winternitz P 2000 Bases for representations of 
quantum algebras J. Phys. A: Math. Gen. 33 5303-13 

[10] Atakishiyev N M and Klimyk A U 2002 Diagonalization of operators 
and one-parameter families of nonstandard bases for representations of 
su g (2) 2002 J. Phys. A: Math. Gen. 35 5267-78 

[11] Gasper G and Rahman M 1990 Basic Hypergeometric Series (Cambridge 
University Press: Cambridge). 

[12] Nomura M 1990 A Jordan-Schwinger representation of quadratic rela- 
tions for SU q (2) operators and the g-analog Wigner-Eckart theorem J. 
Phys. Soc. Japan 59 2345-54; 



19 



[13] Nomura M 1990 An alternative description of the quantum group SU q (2) 
and the g-analog Wigner-Eckart theorem J. Phys. Soc. Japan 59 439- 
48. 

[14] Ballesteros A and Chumakov S M 1999 On the spectrum of a Hamilto- 
nian defined on su g (2) and quantum optical models J. Phys. A: Math. 
Gen 32 6261-9 

[15] N.M. Atakishiyev, S.M. Chumakov and K.B. Wolf 1998, Wigner distri- 
bution function for finite systems, J. Math. Phys. 39 6247-61 

[16] Krawtchouk M 1929 Sur une generalization des polinomes d'Hermite C. 
R. Acad. Set. Pans 189 620-2 

[17] Atakishiyev N M, Pogosyan G S, Vicent L E and Wolf K B 2001 Finite 
two-dimensional oscillator: I. The Cartesian model J. Phys. A: Math. 
Gen. 34 9381-98 

[18] Van der Jeugt J and Jagannathan R 1998 Realizations of su(l, 1) and 
C/ 9 (su(l, 1)) and generating functions for orthogonal polynomials J. 
Math. Phys. 39 5062-78 

[19] Simon R and Wolf K B 2000, Fractional Fourier transforms in two di- 
mensions, J. Opt. Soc. Am. A 17 2368-81. 

[20] Atakishiyev N M, Pogosyan G S, Vicent L E and Wolf K B 2001 Finite 
two-dimensional oscillator: II. The Radial model, J. Phys. A: Math. 
Gen. 34 9399-415 

[21] Atakishiyev N M, Pogosyan G S and Wolf K B 2003 Contraction of the 
finite radial oscillator Int. J. Mod. Phys. A 18 329-41 

[22] Atakishiyev N M, Pogosyan G. S. and Wolf K B 2003 Contraction of 
the finite one-dimensional oscillator Int. J. Mod. Phys. A 18 317-27 

[23] Macfarlane A J 1989 On g-analogues of the quantum harmonic oscillator 
and the quantum group SU(2) q J. Phys. A: Math. Gen. 22 4581-8 

[24] Biedenharn L C 1989 The quantum group SU q (2) and a g-analogue of 
the boson operators J. Phys. A: Math. Gen. 22 L873-L878 



20 



